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I.  INTRODUCTION 


The  ability  to  predict  the  aerodynamics  of  shell  through  the  use  of 
advanced  computational  methods  Is  becoming  an  Important  factor  In  the  design 
an<r  testing  of  current  and  future  systems  of  Interest  to  the  U.S.  Arny. 
State-of-the-art  technology  Is  being  applied  at  BRL  In  a  continuing  research 
effort  directed  towards  computing  the  aerodynamic  forces  that  determine  the 
behavior  of  shell.  Through  advanced  finite-difference  techniques  which  are 
being  used  to  solve  the  tnln-layer  Navi er- Stokes  equations,  improvements  in 
computational  predictive  capabilities  for  shell  are  being  achieved. 

This  paper  documents  a  computational  aerodynamics  stucfy  of  an  Initial 
design  configuration  which  is  Intended  for  use  as  a  training  round  for  a  25mm 
gun.  The  design  configuration  of  interest  is  similar  to  a  standard  spinning 
shell  with  one  important  difference:  the  boattail  is  beveled  such  that 
axl symmetry  no  longer  exists  about  the  projectile  axis.  Referred  to  as  a  non- 
conical  boattail1"2  (NCBT),  this  geometric  feature  causes  a  large  de-spin 
moment  to  act  upon  the  spinning  shell. 

The  non-coni  cal  boattail  is  expected  to  reduce  the  gyroscopic  stability 
enough  to  cause  the  shell  to  become  unstable  in  flight.  Thus  it  is  desirable 
to  determine  the  roll  damping  of  the  shell  as  a  function  of  both  Mach  number 
and  boattail  length.  The  overall  predicted  static  aerodynamic  behavior  of 
this  non-axi symmetric  shell  is  presented.  Comparisons  are  made  with  computa¬ 
tional  predictions ^pr  axisymmetric  configurations  as  well  as  range  firings. 

The  computational  approach  used  in  this  study  is  the  thin-layer  Parabo¬ 
lized  Navier-Stokes  (PNS)  technique  reported  by  Schiff  and  Steger. 3  This  PNS 
technique  has  been  used  successfully  to  predict  the  steady,  t iree-dimensional , 
viscous,  supersonic  flow  for  axisymmetric  spinning  shell  at  small  angles  of 
attack  (a  <  6°). 4-6  This  method  has  also  been  applied  to  non-axi symmetric 
configurations  at  zero  spin  with  encouraging  results.7 

The  PNS  approach  offers  certain  advantages  over  time-dependent  solutions 
that  are  considered  important  for  this  study.  Computer  storage  and  CPU  re¬ 
quirements  are  much  less,  making  parametric  studies  over  a  range  of  Mach  num¬ 
bers  and  angles  of  attack  more  feasible.  At  the  same  time,  it  is  of  interest 
to  apply  the  PNS  technique  to  current  problems  that  test  its  capabilities  and 
define  its  limitations.  This  report  describes  computational  results  using  the 
PNS  technique  for  a  configuration  of  non-axi symmetry  applicable  to  spin  rates 
typical  of  army  gun  twist. 


II.  MODEL  AND  TEST  CONDITIONS 

The  shell  configuration  of  interest  is  depicted  in  Figure  1.  The  shell 
consists  of  a  3  caliber  secant-ogive  nose  and  a  .75  caliber  cylindrical  mid¬ 
section.  The  reference  diameter  of  the  body  is  13.5  mm  (sub-caliber  round), 
occurring  at  the  cylindrical  midsection.  The  afterbody  consists  of  a  triangu¬ 
lar  non-conical  boattail  of  lengths  varying  to  a  maximum  of  2  calibers.  The 
non-conical  boattail  is  made  up  of  three  non-canted  7°  beveled  cuts  spaced  at 
120°  intervals  around  the  body.  The  initial  design  configuration  has  a  1/2 
caliber  non-conical  boattail. 
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The  computational  shell  geometry  possesses  a  sharp  nosetip,  whereas  the 
actual  shell  geometry  possesses  a  rounded  nosetip  with  a  bluntness  ratio  of 
85.  Reference  5  indicates  that  such  a  small  bluntness  ratio  has  a  negligible 
effect  on  the  pitch-plane  aerodynamics  of  shell.  The  blunt  nosetip  contribu¬ 
tion  to  aerodynamic  drag,  however,  is  likely  to  be  noticeable.  In  this  study, 
differences  In  aerodynamic  performance  due  to  contributions  from  other  than 
geometric  non-axi symmetry  are  of  minimal  Interest.  For  purposes  of  clarity, 
however.  It  should  be  noted  that  the  sharp  ogive  nosetip  is  itself  modeled  by 
a  small  conical  extension  possessing  a  half-angle  of  13.7°. 

The  dominant  geometric  feature  of  the  shell  of  Interest  Is  the  triangular 
non-coni  cal  boattail.  Figure  1  shows  the  non-coni  cal  boattail  consisting  of 
three  non-canted  7°  beveled  cuts  spaced  at  120°  intervals  around  the  body. 
Computations  were  also  made  for  two  axisymmetric  configurations  consisting  of 
the  same  forebody  as  in  Figure  1  but  with  0°  and  7°  boattails.  Computations 
were  performed  for  boattail  lengths  varying  from  zero  to  2  calibers.  Figure  1 
shows  a  full  2  caliber  boattail;  shorter  boattail  configurations  were  trunca¬ 
tions  of  the  full  2  caliber  boattail.  Figure  2  shows  the  convention  for  roll 
orientation  for  the  non-axi symmetric  shell. 

The  test  conditions  are  sea-level  atmospheric  .  The  Mach  numbers  are  1.5, 
3,  and  5;  the  Reynolds  numbers  based  on  body  diameter  are  .458  x  106, 
.917  x  lo6  and  1.52  x  io6  respectively.  Both  the  free  stream  static  tempera¬ 
ture  and  the  body  wall  temperature  are  taken  as  273K,  and  the  boundary  layer 
is  assumed  to  be  fully  turbulent.  The  angle  of  attack  is  assumed  to  be  in  the 
linear  range,  from  zero  to  2°.  For  cases  with  spin,  a  non-dimensional  spin 
rate  of  PD/V  *  .1417  corresponds  to  a  gun  twist  of  23.94  calibers  of  transit 
per  revolution  for  this  sub-caliber  -ound. 


III.  COWUTATIONAL  APPROACH 
1.  GOVERNING  EQUATIONS/NUMERICAL  METHOD 

The  results  for  this  study  were  obtained  using  the  thin-layer  Parabolized 
Navier-Stokes  (PNS)  technique  as  formulated  by  Schiff  and  Steger. 3  The  PNS 
equations  can  be  marched  spatially  along  the  body  in  the  direction  where  the 
flow  is  supersonic  and  contains  no  streamwise  separation.  The  thin-layer  ap¬ 
proximation8  retains  the  viscous  derivatives  in  the  radial  direction  only. 

By  employing  the  subsonic  sublayer  approximation  as  discussed  in  Refer¬ 
ence  3,  the  steady,  three-dimensional,  thin-layer  equations  in  strong  conser¬ 
vative  form  can  be  written  as: 

AAA  A 

3E  3F  3G  1  3S 

— L  +  _  +  _  =  ^_  (1) 

3£  3n  3;  Re  3; 


where  S,  n,  c  are  the  fixed  generalized  coordinate  variables  shown  in  Figure 
3,  where: 


C  *  £(x)  Is  the  st reamwlse,  longitudinal  (marching)  coordinate 

n  »  n(x,y,z)  is  the  circumferential  coordinate 

C  *  c(x,y,z)  Is  the  radial  coordinate,  with  the  body 
mapped  onto  the  ;  =  0  plane. 

Equation  (1)  represents  the  thin-layer  approximation  to  the  equations  of 
mass,  momentum,  and  energy  conservation  in  the  three  coordinate  directions. 

^  A  A  A 

The  Invlscld  flux  vectors  Es,  F,  G  and  the  thin-layer  model  viscous  term  S  are 
functions  of  the  dependent  variables  represented  by  the  vector  q»  given  by: 

qT  =  (p,  pu,  pv,  pw,  e) 


Equation  (1)  is  solved  by  applying  the  conservative,  approximately  fac¬ 
tored,  Implicit,  delta  form,  finite-difference  numerical  algorithm  formulated 
by  Beam  and  Warming.9  A  fully  turbulent  boundary  layer  is  simulated  using  the 
two-layer  algebraic  model  of  Baldwin  and  Lomax. 8 

2.  INITIAL  AND  BOUNDARY  CONDITIONS 

An  initial  cross-sectional  plane  of  data  is  obtained  near  the  nosetip  of 
the  projectile  through  the  use  of  a  step-back  procedure  which  employs  the  con¬ 
ical  flow  assumption,  as  outlined  in  References  3,4,  and  10.  Once  the  initial 
plane  of  data  has  been  generated,  the  marching  procedure  can  then  be  applied 
over  the  remainder  of  the  projectile. 

The  outer  boundary,  which  consists  of  the  bow  shock  emanating  from  the 
sharp  nosetip  of  the  projectile,  is  numerically  fit  using  the  implicit  boun¬ 
dary  procedure  of  Reference  10. 

At  the  body  surface,  the  pressure  p$  is  defined  from  =  0  across  the 

subsonic  sublayer,  and  a  no  slip  boundary  condition  is  imposed.  For  the  case 
of  a  projectile  possessing  a  finite  spin  rate  w,  the  velocities  at  the  surface 
become: 

uw  *  0  v„  =  rw  cos  <p  ww  =  -ru>  sin  +  (2) 


where  $  is  the  roll  angle  measured  from  the  windward  side,  using  the  sign 
convention  as  shown  in  Figure  3. 

3.  QUASI-STEADY  APPROXIMATION 

The  governing  equations  for  the  PNS  technique  are  derived  on  the  basis  of 
several  assumptions  about  the  flow  field  (e.g.  supersonic,  thin-layer  approxi¬ 
mation,  et  al).  Of  particular  importance  here  is  the  fact  that  the  flow 
fields  of  References  3  thru  7  share  one  characteristic  which  further  ensures 
the  suitability  of  the  PNS  approach;  that  is,  they  are  steady  flows. 


For  the  present  case  of  flow  over  a  spinning  non-axi symmetric  projectile 
using  a  fixed  coordinate  system,  the  steady  nature  of  the  flow  is  lost.  Even 
at  o  *  0°,  a  constant  change  in  roll  orientation  of  the  projectile  contributes 
a  periodic  unsteadiness  to  the  flow  field.  In  effect,  this  is  an  unsteady 
problem  being  modeled  using  a  steady  approach.  Particular  attention  must  be 
focused  on  the  unsteady  nature  of  this  flow  and  the  ramifications  of  its 
modeling  using  a  steady  approach. 

To  apply  the  PNS  technique  to  the  case  of  a  spinning  non-axi  symmetric 
bo<fy  is  to  assume  that  the  flow  is  quasi-steady.  Such  an  assumption  Implies 
that  for  some  instant  of  time,  the  rates  of  change  of  the  properties  at  every 
point  in  the  computational  domain  are  nearly  zero.  This  is  a  reasonable 
assumption  if  the  nondimensional  spin  rate,  PD/V,  is  small  enough. 

A  series  of  computations  were  run  to  determine  a  range  of  nondimensional 
spin  rates  over  which  the  quasi-steady  assumption  was  valid.  Cases  were  run 
varying  PD/V  from  small  values  to  typical  flight  values.  As  shown  in  Figure 
5,  the  computed  roll  damping  was  relatively  insensitive  to  spin  rate  for  the 
range  .005  <  PD/V  <  .05.  The  linear  behavior  of  rolling  moment  with  respect 
to  spin  rate  allows  the  computed  roll  damping  coefficient  to  be  applied  to  the 
much  higher  spin  rates  that  are  typical  of  launch  and  flight  conditions  for 
shell. 


IV.  PRELIMINARY  ANALYSIS 


1.  BACKGROUND 

The  current  objective  in  the  design  of  a  training  round  with  a  non- 
conical  boattail  is  to  reduce  its  firing  range.  The  larger  de-spin  moment 
acquired  with  a  non-conical  boattail  is  expected  to  induce  gyroscopic 
instability  faster  than  would  a  completely  axisymmetric  configuration. 
Gyroscopic  instability  is  associated  with  a  fast-growing  average  amplitude  of 
angle  of  yaw. 

For  a  symmetric  body  of  revolution,  the  total  drag  behaves  as  a  linear 
function  of  the  square  of  the  angle  of  yaw. 11  The  general  effect  of  the  non- 
tonical  boattail  on  the  spinning  shell  should  therefore  be  increased  drag  due 
to  the  larger  angles  of  attack  experienced  during  flight.  It  is  of  primary  in¬ 
terest  in  this  study  to  determine  the  characteristics  of  such  a  configuration 
at  small  angle  of  attack. 

Some  aerodynamic  effects  of  shell  with  non-axi symmetric  boattails  have 
previously  been  examined  and  give  useful  information  in  terms  of  what  kind  of 
performance  to  expect  from  a  triangular  non-conical  boattail.  Platou12"13 
examined  an  array  of  projectiles  with  non-axi symmetric  boattails  in  order  to 
determine  the  merits  of  such  shapes.  Wind  tunnel  tests  and  range  firings  of 
configurations  with  non-conical  boattails  consistently  indicated  improved 
static  stability  over  configurations  with  conical  boattails. 

The  flat  faces  of  the  straight  non-conical  boattails  were  described  in 
Reference  12  as  "fin-like  surfaces"  possessing  no  differential  cant.  As  such, 
boattails  similar  to  the  design  of  Interest  in  this  study  were  shown  to  have 
much  higher  roll  damping  than  their  axisymmetric  counterparts.  A  two-caliber, 
7°  triangular  boattail  on  a  5  to  7  caliber  shell  was  shown  through  free-flight 

4 


testing  at  transonic  speeds  to  have  a  roll  damping  coefficient  an  order  of 
magnitude  larger  than  a  cylindrical  boattall.* 

Reference  12  also  showed  a  significant  dependence  of  the  rolling  moment 
on  roll  orientation  and  angle  of  attack.  For  the  same  configuration  just  men¬ 
tioned  C£  was  found  to  be  of  the  same  order  of  magnitude  as  C£  . 

°  P 

More  recently,  wind  tunnel  roll -oscillation  tests14  were  performed  on  a 
shape  quite  similar  to  the  one  of  Interest  in  this  study.  The  primary  objec¬ 
tive  of  the  tests  was  to  determine  the  dynamic  stability  characteristics  in 
roll.  The  total  length  of  the  configuration  was  6  calibers  and  the  length  of 
the  7°  triangular  non-coni  cal  boattall  was  1  caliber.  The  large  degree  of  un¬ 
certainty  present  in  the  reduced  data  (typically  50%  of  the  total  magnitude) 
makes  the  roll  characteristics  difficult  to  Interpret  qualitatively.  At  Mach 
2,  the  roll  damping  was  practically  independent  of  boattail  roll  orientation; 
likewise,  variation  of  roll  damping  did  not  occur  as  a  function  of  angle  of 
attack  for  a  <  10°.  The  pitch-plane  aerodynamics  also  showed  little  or  no 
dependence  on  rcll  orientation,  and  in  this  sense  the  model  was  viewed  as 
being  nearly  ax i symmetric. 


2.  STATIC  AERODYNAMIC  ANALYSIS 

Sign  conventions  for  all  aerodynamic  coefficients  discussed  here  are 
based  on  the  coordinate  system  of  Reference  11,  shown  in  Figure  4. 

According  to  Murphy11,  the  rolling  moment  coefficient  for  a  non-axisym- 
metric  configuration  can  be  written  as: 


C£  -if  C  +  (PD/V)  C  +  6C  .  (3) 

6  p  a 


It  can  be  seen  from  Figure  2  that  the  triangular  non-coni  cal  boattail  is  ana¬ 
logous  to  a  set  of  three  control  surfaces  with  no  differential  cant;  there¬ 
fore  C.  is  zero.  It  remains  then  to  determine  the  values  of  C  and  C  for 

X  •.  X  X 

o  pa 

the  flight  conditions  of  interest.  If  the  yaw  angle  is  zero,  then  the  rolling 

moment  coefficient  can  be  described  using  the  roll  damping,  C,  ,  alone. 

P 

Similarly,  if  the  projectile  is  not  spinning,  then  the  rolling  moment  coeffi¬ 
cient  can  be  described  using  C,  alone.  If  the  projectile  is  spinning  at 

a 

some  yaw  angle,  both  remaining  terms  on  the  right-hand  side  of  Equation  (3) 
must  be  considered  in  evaluating  the  rolling  moment. 


*  It  is  of  interest  to  note  that  a  major  design  objective  in  References  12  and 
13  was  to  increase  the  range  and  accuracy  of  shell  by  employing  a  non- 
conical  boattail.  In  order  to  maintain  the  gyroscopic  stability  of  that 
family  of  shell,  the  next  design  change  was  to  cant  or  twist  the  boattail  to 
correspond  to  the  gun  barrel  twist. 
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The  simplest  assumed  form  of  C,  for  an  n-finned  missile  that  satisfies 
the  Mapl e-Synge 15  requirement  of  analytic  continuity  in  both  pitch  and  roll  is 

C-  ■  A  sin  n6 

(4) 

where  A  =  k6n”* ,  k  *  constant. 

For  an  axisymmetric  body  of  revolution,  CA  =  0.  For  the  projectile  of 

a 

interest  in  this  study,  n  =  3;  therefore  the  rolling  moment  contribution  due 
to  non-axi symmetry  alone  is  periodic  in  120°  roll  intervals.  From  Figure  2  and 
Equation  (4),  it  is  apparent  that  C,  is  known  for  all  roll  orientations  if  it 
is  known  for  0°  <  e  <  60°.  a 

More  importantly,  if  such  an  assumed  form  for  is  valid,  then  at  a 

a 

constant  spin  rate  must  oscillate  about  zero  at  a  frequency  which  is  pro¬ 
b'd 

portional  to  the  spin  rate.  (We  note  that  the  PNS  computation  is  carried  out 

with  a  zero  coning  rate,  therefore  P  =  8.)  If  the  spin  rate  is  large  enough, 
then  the  roll  frequency  becomes  large  compared  to  the  pitch  and  yaw  frequen¬ 
cies,  and  the  integrated  effect  of  C£  along  the  trajectory  becomes  nil. 


The  spin  rates  of  spin-stabilized  shell  are  normally  quite  large  compared 
to  those  of  fin-stabilized  projectiles.  The  determination  of  a  minimum  spin 
rate  at  which  the  effect  of  C.  can  no  longer  be  ignored  in  the  prediction  of 

*o 

projectile  flight  is  unclear;  the  judgement  concerning  the  importance  of  the 
effect  at  various  spin  rates  is  not  addressed  here.  It  is,  however,  desirable 
to  establish  the  relative  magnitudes  of  such  of  oscillations  compared  to  the 
nominal  values  of  the  computed  aerodynamic  coefficients.  Section  V  will  dis¬ 
cuss  further  the  effect  of  roll  orientation  on  the  static  aerodynamic 
coefficients. 

The  PNS  code  is  employed  to  compute  the  flow  field  over  the  shell  by 
spacially  marching  the  solution  downstream.  In  general,  the  solution  is  no 
longer  valid  when  regions  of  axial  separation  or  subsonic  flow  are  en¬ 
countered.  Therefore,  the  domain  of  the  PNS  computation  does  not  include  the 
base  region.  The  rolling  and  pitch-plane  aerodynamic  characteristics  of  shell 
are  not  overly  sensitive  to  base  effects  in  the  supersonic  flight  regime.  The 
same  cannot  be  said  about  axial  force  since  base  drag  at  the  Mach  numbers  of 
Interest  is  a  significant  component  of  the  total  drag. 

For  these  reasons  a  second  computer  code  -  MC  DRAG16  -  is  used  to  aid  in 
the  analysis  of  drag  for  the  shell  of  interest.  MC  DRAG  is  a  fast  running, 
easy  to  use  program  which  employs  empirical  data,  low  order  approximations, 
and  similarity  concepts  to  estimate  the  components  of  zero-yaw  drag  as  fun¬ 
ctions  of  Mach  number.  For  this  study,  it  is  of  main  interest  to  apply 
MC  DRAG  to  obtain  estimates  of  base  drag  for  the  non-axi symmetric  shapes. 


Following  the  notation  of  Reference  16,  the  zero-yaw  drag  coefficient 
takes  the  form: 


C°0  *  C°H  *  St  *  \  *  V 


Furthermore,  the  total  forebody  drag  coefficient  is  defined  as: 


^D  =  ^D  +  ^D  +  ^n  • 
UF  UH  UBT  USF 


In  order  to  obtain  base  drag  estimates  for  the  non-axi symmetric  shapes 
from  MC  DRAG,  a  simple  geometric  model  for  the  non-conical  boattail  must  be 
established.  This  is  necessary  because  MC  DRAG  can  only  be  applied  to 
boattails  which  are  conical. 

Figure  6  shows  the  cross-sectional  areas  as  functions  of  axial  position 
for  two  of  the  boattails  under  consideration.  Each  configuration  consists  of 
the  same  nose  and  cylindrical  midsection,  differing  only  in  the  boattail  geo¬ 
metry.  One  is  a  7°  conical  boattail  and  the  other  is  the  triangular  non- 
conical  boattail  described  in  Section  II. 


It  is  well  known  that 
sectional  area  of  the  base. 


boattail  should  be  bounded  by  Cn  for  the  7°  concial  boattail  and  the 

UB 

cylindrical  extension.  Thus,  the  base  drag  for  the  non-axi symmetric 
configuration  has  been  estimated  from  the  drag  of  a  conical  boattail  which 
possesses  an  equivalent  base  area.**  No  further  distinction  is  made  between 
the  non-axi symmetric  shape  and  this  simple  geometric  model  in  the  ensuing 
discussion  of  drag. 


bounded 


base  drag  is  a  strong  function  of  the  cross- 
For  a  given  length,  tn  for  the  non-conical 

B 

by  Cn  for  the  7°  concial  boattail  and  the 


V.  RESULTS 

Solutions  were  generated  over  both  conical  and  non-conical  boattails 
using  an  algebraic  grid.  The  grid  consisted  of  45  exponentially-stretched 
points  in  the  radial  direction  from  the  body  to  the  shock  and  spacing  of 
points  circumferentially  around  the  body  at  5°  intervals.  Grid  spacing  at  the 
wall  was  maintained  such  that  the  first  point  above+the  body  was  contained  in 
the  laminar  sublayer  (specifically,  a  value  of  3  <  y  <5  was  used  to  maintain 
adequate  viscous  resolution). 


**  Reference  1  examined  an  "equivalent  area"  boattail  shape  to  model  the  non- 
conical  boattail,  unlike  the  geometric  model  employed  here.  The  equivalent 
area  boattail  was  an  axisymmetric  configuration  whose  cross-sectional  area 
was  equal  to  that  of  the  non-conical  boattail  at  all  axial  locations. 
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All  numerical  computations  were  performed  on  a  CDC  7600  computer  with  a 
speed  of  .0026  CPU  sec/step/grid  point  (number  of  grid  points/step  =  45  x 
72).  The  computer  time  required  to  obtain  a  starting  solution  and  march  to  an 
axial  position  just  upstream  of  the  boattail  was  typically  2000  CPU  seconds. 
Another  1000  CPU  seconds  were  required  to  march  over  each  of  the  3  different 
boattail  configurations. 

1.  ROLL  DAMPING 

Figure  7  shows  the  roll  damping  coefficient,  C*  ,  developing  over  the 

P 

shell  of  Interest  at  Mach  5.  Comparison  is  made  between  the  non-coni  cal 
boattail  configuration  and  the  two  axi symmetric  configurations.  It  is  apparent 
that  the  full  2  caliber  non-coni  cal  boattail  generates  an  order  of  magnitude 
increase  in  roll  damping  compared  to  the  axi symmetric  configurations.  The 
order  of  magnitude  increase  In  roll  damping  for  the  non-conical  configuration 
is  almost  entirely  attributable  to  surface  pressure  rather  than  viscous 
effects. 

Figure  8  shows  the  development  of  C.  over  the  non-axi symmetric  shell  for 

P 

Mach  numbers  1.5,  3  and  5.  It  is  apparent  that  the  effect  of  C.  decreases 

P 

with  increasing  Mach  number  in  this  speed  regime.  Slight  oscillations  in  the 
predicted  roll  damping  occur  for  the  non-axi symmetric  bocty.  These  oscilla¬ 
tions  are  a  result  of  grid  coarseness  in  the  circumferential  direction,  which 
affects  the  solution  near  the  parabolic  intersections  of  the  boattail  cone  and 
beveled  flats.  Figure  9  presents  the  roll  damping  as  a  function  of  Mach 
number  using  boattail  length  as  a  parameter. 

Figure  10  compares  the  computed  C.  to  that  obtained  from  range  firings17 

P 

for  the  1/2  caliber  non-conical  boattail  configuration.  Range  reductions 
showed  that  not  all  experimental  data  points  were  wel 1 -determined,  therefore 
some  scatter  exists.  Disagreement  between  computation  and  experiment  below 
Mach  3  is  from  0 %  to  20%,  while  disagreement  near  Mach  5  is  about  30%. 

2.  PITCH-PLANE  AERODYNAMICS 

Figure  11  shows  the  development  of  normal  force  coefficient  over  the 
shell  for  conical  and  non-conical  boattails  at  Mach  1.5  .  Only  slight  differ¬ 
ences  exist  in  the  normal  force  due  to  the  three  different  boattails.  Figure 
12  shows  the  small-a  slope  of  the  normal  force  coefficient  for  the  non-conical 
configuration  as  a  function  of  Mach  number  using  boattail  length  as  a  para¬ 
meter.  Figure  13  uses  the  same  approach  to  present  the  pitch-plane  center  of 
pressure.  Both  CN  and  CP  are  seen  to  be  non-linear  with  respect  to  boattail 
length.  a 

Pitching  moment  coefficients  are  presented  here  with  reference  to  the 
centers  of  gravity  of  each  individual  configuration.  Table  1  shows  the  c.g. 's 
for  the  various  configurations,  assuming  solid,  isotropic  material.  The  only 
addition  is  a  .5  caliber  non-conical  boattail  case,  possessing  a  small  base 
cavity,  and  used  in  the  firing  range  comparison. 
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Figure  14  shows  the  small -a  slope  of  the  pitching  moment  coefficient  for 
the  non-axi symmetric  configuration  as  a  function  of  Mach  number,  using  boat- 
tail  length  as  a  parameter.  It  is  apparent  that  Cm  decreases  with  increasing 

a 

Mach  number,  and  increases  somewhat  linearly  with  increasing  boattail  length. 

Figure  15  compares  the  computed  Cm  to  that  obtained  from  range  firings 

"a 

for  the  1/2  caliber  non-conical  boattail  configuration  with  small  base  cavity. 
Disagreement  between  computation  and  experiment  is  at  most  5%. 

Figures  16,  17,  18  and  19  compare  values  of  pitching  moment  coefficient 
for  configurations  having  conical  and  non-conical  boattails,  keeping  boattail 
length  fixed.  For  a  boattail  of  length  1/2  caliber,  differences  in  pitching 
moment  due  to  boattail  shape  can  barely  be  discerned.  As  boattail  length  in¬ 
creases,  these  differences  become  more  apparent.  In  general,  the  fin-like 
Influence  of  the  non-conical  boattail  is  demonstrated  by  the  more  stable 
values  of  pitching  moment  compared  to  the  axisymmetric  shapes. 

3.  DRAG 


Figures  20  and  21  show  the  total  forebody  drag  developing  over  the 
axisymmetric  and  non-axi symmetric  shell  for  Mach  numbers  1.5  and  3,  respec¬ 
tively.  Comparison  is  made  with  the  predicted  values  from  MC  DRAG.  For  each 

case,  Cn  for  the  non-axi symmetric  shape  falls  between  Cn  for  the  two 
uF  uF 


axisymmetric  shapes. 


Figure  22  shows  Cn  from  the  PNS  code  for  the  non 
UF 


axisymmetric  shape  as  a  function  of  Mach  number  using  boattail  length  as  a 
parameter.  Figure  23  shows  the  base  drag  for  the  non-axi  symmetric  shape 
obtained  from  MC  DRAG  as  a  function  of  Mach  number,  using  boattail  length  as  a 
parameter. 


Figure  24  shows  the  predicted  total  drag  of  the  non-axi symmetric  shell  as 
a  function  of  Mach  number,  using  boattail  length  as  a  parameter.  The  total 
drag  is  obtained  by  combining  the  forebody  drag  from  PNS  and  the  base  drag 
from  MC  DRAG.  It  is  apparent  that  the  total  drag  decreases  with  increasing 
Mach  number  and  decreases  slightly  with  increasing  boattail  length. 


Figure  25  compares  the  predicted  total  drag  with  that  obtained  from  range 
firings.  Disagreement  below  Mach  3  is  within  10%,  while  disagreement  near 
Mach  5  is  about  20%. 


4.  ROLL  ORIENTATION 

The  effect  of  roll  orientation  on  the  roll  and  pitch-plane  aerodynamic 

coefficients  is  summarized  in  Table  2.  Cases  were  run  for  5  roll  orientations 

at  2°  angle  of  attack  and  no  spin.  For  non-conical  boattail  lengths  of  1  and 

2  calibers,  the  computations  showed  little  effect  of  roll  orientation  on 

C,,  ,  CP,  or  C  .  The  largest  variation  of  Cm  was  about  5%,  occurring  for 
N  m  3  m 

a  a  a 

the  case  of  Mach  1.5  with  a  2  caliber  non-conical  boattail. 
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Table  2  also  shows  the  computed  rolling  moments  for  each  case.  Using 
these  results  and  the  form  for  already  discussed  in  Section  IV,  the  ampli- 

*a 

tude  A  of  the  sinusoidal  rolling  moment  was  obtained  for  each  case  and  is 
shown.  Note  that  the  form  for  employed  here  takes  roll  orientation  into 

*a 

account;  therefore,  the  amplitude  A  of  at  a  given  Mach  number,  angle  of 

a 

attack,  and  boattail  length  should  be  equal  for  all  roll  orientations.  From 
symmetry,  the  rolling  moments  for  6=0°  and  8  =  60°  should  equal  zero,  and 
the  maximum  rolling  moment  should  occur  at  0  *  30°.  It  is  apparent  that  the 
values  for  and  A  in  Table  2  do  not  satisfy  Equation  5  adequately.  Two 

possible  reasons  are  (1)  Equation  (4)  is  an  over-simplified  form  for  Ct  ,  or 

a 

(2)  numerical  damping  in  the  PNS  solution  Inhibits  precise  determination 
of  C.  for  a  nearly  axisymmetric  body, 
a 

From  the  results  presented  thus  far  and  the  form  for  discussed  in  Sec¬ 
tion  IV,  it  can  be  deduced  that  even  at  moderate  angles  of  attack  the  rolling 
moment  induced  by  on  the  flight  of  the  projectile  is  negligible  compared 

a 

to  that  which  is  induced  by  C  .  At  Mach  1.5,  with  a  1  caliber  non-conical 

P 

boattail  and  a  non-dimensional  spin  rate  PD/V  =  .1417  (launch  condition),  the 
computed  t  gives  a  rolling  moment  of  approximately  C  *  -.0085.  At  a=  10° 
P 

the  computed  (disregarding  non-linear  effects)  gives  a  maximum  rolling 
a 

moment  of  approximately  =  -.00021,  at  least  an  order  of  magnitude  less  than 

the  contribution  due  to  . 

P 

Similarly,  the  maximum  side  moment  produced  by  geometric  non-symmetry 
about  the  pitch-plane  was  compared  to  Magnus  moment  at  each  Mach  number. 
Following  the  same  form  as  Equation  5,  the  computed  maximum  geometric  side 
moment  occured  for  e  =  30°.  For  the  1/2  caliber  non-conical  boattail 
configuration,  the  maximum  geometric  side  moment  was  computed  to  be  at  least 
an  order  of  magnitude  less  than  the  Magnus  moment  at  launch  spin  rate. 

It  was  concluded  that  roll  orientation  has  a  negligibly  small  effect  on 
pitch  plane,  roll  damping  and  Magnus  coefficients  for  the  cases  where  the  NCBT 
is  less  than  one  caliber.  Variations  of  these  coefficients  with  respect  to 
roll  orientation  were  Ignored  In  the  trajectory  simulations  which  are 
presented. 

5.  TRAJECTORY  SIMULATION 

In  order  to  aid  in  the  design  of  a  range-limited  projectile,  the  computed 
aerodynamic  coefficients  were  applied  to  a  mathematical  simulation  of  the 
projectile  in  flight.  A  six  degree-of-freedom  computer  code  was  used  to 
Integrate  the  equations  of  motion  for  a  rigid  projectile.18  Using  as  input 


the  aerodynamic  forces  and  moments  computed  In  this  study,  ranges  at  which  the 
gun-launched  non-ax 1  symmetric  configurations  would  experience  large  Increases 
In  yaw  angle  were  determined. 

The  trajectory  simulation  code  requires  the  physical  properties  of  the 
projectile  (see  Table  1)  and  Initial  flight  conditions.  In  addition,  an  array 
of  aerodynamic  coefficients  must  be  specified  as  functions  of  Mach  number  and 
angle  of  attack.  Some  coefficients,  such  as  pitch  damping  (Cm  +  Cm.) 

d  o 

and  cubic  pitching  moment  (Cm  ),  may  have  a  significant  effect  on  the  trajec- 

a3 

tory,  but  are  not  directly  available  and  must  be  estimated.  The  overall  ef¬ 
fect  of  such  estimations  was  gauged  through  parametric  variations  with  guid¬ 
ance  from  Reference  19. 

The  launch  velocity  Is  taken  to  be  Mach  4. 78. 20  The  gun  twist  Is  taken 
to  be  23.94  calibers/revolution20  and  the  launch  elevation  Is  set  at  35.55 
mils  (2°).  The  initial  angle  of  attack  Is  taken  to  be  zero.  The  Initial 
pitch  rate  is  taken  to  be  10  radians/second.  In  order  to  provide  a  first 
maximum  yaw  between  1°  and  2°. 

Results  from  trajectory  simulations  are  presented  here  for  two 
configurations:  1/2  caliber  and  3/4  caliber  non-conical  boattail  shell.  As 
previously  mentioned,  the  total  lengths  of  these  two  shell  differ  by  1/4 
caliber. 

Figure  26  shows  Mach  number  as  a  function  of  range  for  the  two  config¬ 
urations.  Both  shell  maintained  forward  flight  long  enough  to  reach  Mach  1, 
occurring  at  ranges  between  2000  and  2500  meters.  Figure  27,  which  shows 
total  yaw  angle  vs.  range,  reveals  very  different  behavior  between  the  two 
shell.  The  1/2  caliber  non-conical  boattail  shell  experienced  small  growth  in 
yaw  over  the  entire  trajectory,  whereas  the  3/4  caliber  non-conical  boattail 
shell  reached  yaw  angles  of  8°  by  1000  meters.  It  is  clear  from  Figures  26 
and  27  that  the  large  yaw  experienced  by  the  longer  shell  caused  significant 
loss  in  velocity  due  to  the  higher  drag. 

Figure  28  shows  nondimensional  spin  rate  PD/V  vs.  range.  The  1/2  caliber 
boattail  shell  experienced  a  continual  increase  in  PD/V  out  to  2500  meters. 
The  3/4  caliber  boattail  shell,  due  to  higher  roll  damping,  maintained  a 
relatively  constant  PD/V. 

Figure  29  shows  gyroscopic  stability11  vs.  range.  At  launch,  both  shell 

possess  Sg  <  1.5,  which  is  a  nominal  value  for  spin-stabilized  projectiles. 

The  1/2  caliber  boattail  shell  experienced  a  continual  increase  in  Sg  out  to 

2500  meters.  The  3/4  caliber  boattail  shell  experienced  a  decrease  in  Sg, 

falling  below  1.0  at  a  range  of  about  600  meters,  and  climbing  above  1.0 
through  the  remainder  of  the  trajectory. 
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VI.  CONCLUSIONS 


A  computational  aerodynamics  study  has  been  made  for  an  initial  design 
configuration  intended  as  a  sub-caliber  training  round  for  a  25mm  gun.  The 
configuration  has  a  triangular  non-conical  boattail  which  Induces  a  large  de¬ 
spin  moment.  The  thin-layer  Parabolized  Navler-Stokes  computational  technique 
was  used  to  compute  the  static  aerodynamic  coefficients  acting  on  the  shell  at 
supersonic  speeds. 

The  computations  showed  that  a  full  2  caliber,  7°,  triangular  non-conical 
boattail  increases  the  roll  damping  by  more  than  an  order  of  magnitude 
compared  to  axi symmetric  configurations  of  the  same  caliber  and  length.  The 
computed  roll  damping  for  the  case  of  a  1/2  caliber  non-conical  boattail  is 
almost  double  that  of  an  axlsymmetric  configuration  of  the  same  caliber  and 
length. 

The  pitch  plane  results  verified  that  the  non-conical  boattail  gives  the 
shell  a  more  statically  stable  behavior  compared  to  axi symmetric  configura¬ 
tions.  The  computed  drag  of  the  non-axi symmetric  shell  was  consistently 
bounded  by  the  drag  of  a  shell  with  a  0°  boattail  and  a  shell  with  a  conical 
boattail  with  the  same  slope  as  the  beveled  cuts  as  the  non-conical  boattail. 
The  effect  of  roll  orientation  on  the  pitch-plane,  Magnus,  and  rolling  forces 
and  moments  was  found  to  be  negligible  for  small  angles  of  attack. 

Comparison  of  computed  aerodynamic  coefficients  to  recently  obtained 
firing  range  data  yielded  the  following:  (1)  roll  damping  coefficient  agreed 
from  0%  to  20%,  except  near  Mach  5,  where  agreement  was  within  30%;  (2)  linear 
pitching  moment  coefficient  agreed  within  5%;  and  (3)  zero  yaw  drag  coeffi¬ 
cient  agreed  with. n  10%,  except  near  Mach  5,  where  agreement  was  within  20%. 

Six  degree-of-f reedom  trajectory  simulations,  based  on  computationally 
determined  static  aerodynamic  coefficients,  were  computed  for  1/2  and  3/4 
caliber  non-conical  boattail  configurations.  The  1/2  caliber  non-conical 
boattail  shell  maintained  small  angles  of  yaw  past  2500  meters  in  range.  The 
3/4  caliber  boattail  shell  reached  moderate  angles  of  yaw  within  a  range  of 
1000  meters.  Since  computed  aerodynamic  coefficients  were  not  available  below 
Mach  1,  no  assessment  was  made  of  the  flight  performance  of  the  shell  beyond  a 
range  of  2500  meters. 


TABLE  1.  Physical  Characteristics  of  Configurations 
of  Interest;  Material:  Steel. 


MOMENTS  OF  INERTIA 
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BOATTAIL 

BOATTAIL 

C.G. 
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AXIAL 

TRANSVERSE 
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CONFIGURATION 

(cal  from  nose) 

(g) 

!x 

!y 

Jt 

3.75 

0 
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2.60 

30.7 

5.53 

37.0 

4.00 

0.25 

7°  non-conlcal 

2.74 

34.5 

6.37 

47.4 

4.25 

0.5 

0° 
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38.3 

7.25 

59.9 

fi 

V 

7°  conical 
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37.4 

6.88 

48.5 

V 

7°  non-conical 

2.87 

38.1 

7.16 

59.2 

V 

\ 

7°  non-conical 

2.78 

_  _ _ 

____ 

—  —  _  _ 

’■ 

w/base  cavity 

i* 

S 

i 

4.50 

0.75 

7°  non-conical 

2.99 

41.5 

7.86 

72.4 

a 

4.75 

1.0 

0° 

3.14 

45.9 

8.98 

91.4 

s 

7°  conical 

3.04 

42.5 

7.65 

64.0 

7°  non-conical 

3.11 

44.6 

8.47 

86.8 

5.25 

1.5 

0° 

3.41 

53.5 

10.71 

133.4 

«; 

7°  conical 

3.20 

46.1 

8.06 

82.0 

i 

7°  non-conical 

3.31 

49.9 

9.41 

118.! 

< 

5.75 

2.0 

0° 

3.67 

61.0 

12.43 

187.! 

% 

7°  conical 

3.31 

48.6 

8.24 

100.2 

it 

7°  non-conical 

3.47 

53.9 

9.97 
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TABLE  2.  Effect  of  Roll  Orientation  on  Various 
Static  Aeroctynamlc  Coefficients. 


A.  1  Caliber  Non-Coni  cal  Boattall 

Case,  6  ■  2° 

,  PD /V  -  0 

Roll 

Orientation 

e 

Mach  No. 

-  1.5 

2° 

ii° 

30° 

ii° 

%  var 

CN 

a 

2.75 

2.76 

2.78 

2.81 

2.82 

2.5% 

CP 

2.21 

2.21 

2.23 

2.25 

2.26 

2.2% 

Cm 

2.46 

2.46 

2.43 

2.40 

2.38 

3.3% 

a 

C* 

3.8*1 0"9 
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— 

A 

— 
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-.00069 

— 

— 

Mach  No. 

=  3 

‘"a 
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0.0% 

CP 

2.37 
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0.0% 

Cm 
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2.10 

2.10 
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3.0*10* 1 1 
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— 

A 

— 
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— 
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CN 
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3.4x1 0*  1 1 

— 

A 
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.00032 
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— 

— 
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TABLE  2.  Effect  of  roll  Orientation  on  Various 
Static  Aerodynamic  Coefficients.  (Continued) 

B.  2  Caliber  Non-Coni  cal  Boattail  Case.  S  •  2°,  PD/V  *  0 
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Roll 
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Mach 
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CN 
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2.71 

2.77 
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Cm 
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CN 
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a 
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Lt 
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A 
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— 

— 
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°N 
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CP 
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a 
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1.1% 

C* 
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— 
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Figure  5.  Validation  of  quasi-steady  assumption  for  roll  damping,  Mach 
a  =  0°,  nnn-axi symmetric  configuration. _ 
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:tional  areas  of  triangular  non-conical  boattail  and  7 
boattail  as  functions  of  axial  position. _ 
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Roll  damping  of  1/2  caliber  non-conlcal  boattail  configuration;  comparison 
between  computation  and  experiment. _ 
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0°  BOATTAIL 

NON-CONICAL  BOATTAIL 


Figure  11.  Normal  force  coefficient  developing  over  body,  Mach 
axi symmetric  and  non-axl symmetric  configurations. 


fferent  boattai  1  lengths. 


Figure  15.  Linear  pitching  moment  coefficient  for  1/2  caliber  non-conical  boattail  configuration 
with  small  base  cavity;  comparison  between  computation  and  experiment. _ 
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moment  coefficient  for  axisymmetric  and  non-axi symmetric 
with  1/2  caliber  boattail. 


pitching  moment  coefficient  for  axi symmetric  and  non-axi symmetric 
urations  with  2  caliber  boattail.  _ _ 


Figure  20.  Forebody  drag  coefficient  developing  over  body,  Mach 
axi symmetric  and  non-axi symmetric  configurations. 
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Forebody  drag  coefficient  for  non-axi symmetric  configuration  with 
different  boattail  lengths,  a=  0°. 
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Figure  25.  Total  drag  coefficient  for  1/2  caliber  non-conical  boattail  configuration; 

comparison  between  computation  and  experiment. _ 
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Trajectory  simulation  for  non-axi symmetric  conf Igu rat  ions;  nondimensional 
spin  rate  as  a  function  of  range. 
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LIST  OF  SYMBOLS 


=  amplitude  of  with  respect  to  roll  orientation 
a 

=  speed  of  sound 

=  total  drag  coefficient  at  a  =  0° 

=  pressure  drag  coefficient  due  to  projectile  base 

=  pressure  drag  coefficient  due  to  projectile  boattail 

*  drag  coefficient  due  to  projectile  forebody 

=  pressure  drag  coefficient  due  to  projectile  head  (nose) 

=  skin  friction  drag  coefficient  due  to  entire  projectile 
wetted  surface 

=  rolling  moment  coefficient 

=  slope  of  rolling  moment  coefficient  due  to  angle  of  attack 
and  roll  orientation 

=  3C£/3(PD/V),  roll  damping  coefficient 


=  3CA/3(6f)  slope  of  roll  moment  coefficient  due  to  canted 

or  deflected  lifting  surfaces 

Cm  =  linear  pitching  moment  coefficient 

a 

Cm  =  cubic  pitching  moment  coefficient 

a  3 

Cm  +  ^m*  =  Pitch  damping  coefficient 

q  a 

CN  =  dC^/da,  zero-a  slope  of  normal  force  coefficient 

a 

CP  =  center  of  pressure  of  normal  force  coefficient 

D  =  diameter  of  model 

e  =  total  energy  per  unit  volume  of  fluid,  normalized  by  p^a.2 
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Es,  F,  G 


u,v,w 


x  f  y  >  z 


LIST  OF  SYMBOLS  (Continued) 

=  flux  vectors  of  transformed  gasdynamic  equations 
=  axial  moment  of  inertia 
=  transverse  moment  of  Inertia 
=  Mach  number 

=  twist  (spin)  of  model  in  calibers  per  revolution 
=  number  of  fins 

=  spin  rate  about  model  axis  in  rad/sec  relative  to  cross-flow 
=  non-dimensional  spin  rate  about  model  axis 
=  pressure  normalized  by 
=  subsonic  sublayer  pressure 
=  vector  of  dependent  variables 
=  Reynold's  number, 

=  local  distance  from  model  axis  to  body  surface 

=  gyroscopic  stability  factor  (unconditionally  unstable  if 
Sg  <  1) 

=  viscous  flux  vector 

=  Cartesian  velocity  components  along  x,y,z  axes,  respectively, 
normalized  by  a*,. 

=  free  stream  velocity 

=  physical  Cartesian  coordinates 

=  distance  from  wall  in  law  of  the  wall  coordinates 

=  angle  of  attack 

=  angle  of  sideslip 

=  angle  of  yaw,  |  ✓  a2  +  02  j 

=  angle  of  fin  deflection  or  cant 


$ 

v 

t 


8 
i 


S,  n,  c 


Subscripts 


Superscripts 


LIST  OF  SYMBOLS  (Continued) 


computational  coordinates  in  the  axial,  circumferential, 
and  radial  directions 


=  roll  orientation  of  shell  relative  to  cross-flow 


=  coefficient  of  viscosity,  normalized  by  free  stream  value 


=  density,  normalized  by  free  stream  value 


roll  angle;  circumferential  position  on  the  body  measured 
from  windward  side 


=  spin  rate  about  model  axis  in  rad/sec 


free  stream  conditions 


sonic  layer 


body  surface  value 


=  transpose 
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This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comments /answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ _ _ Date  of  Report _ 

2.  Date  Report  Received _ _ 

3.  Does  this  report  satisfy  a  need?  (Coament  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Inforaation  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ ' _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments,  What  do  you  think  should  be  changed  to  iaprove  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


